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Remotely  sensed  bathymetry  from  multispectral  imagery 
generally  is  based  on  a  simple  reflectance  model1;  the  radi¬ 
ance  in  wavelength  band  i  at  water  depth  Z  is  given  by  the 
equation 

L,  =  +  c,  •  Rm  •  exp(-2A,  •  Z), 

where  L,  is  the  radiance  value  in  band  i,  Li .  is  the  average 
signal  over  deep  water,  c,  is  a  constant  that  is  a  function  of 
several  optical  parameters,  Ra,  is  the  bottom  reflectance  in 
band  i  over  bottom  type  a,  and  k,  is  the  diffuse  attenuation 
coefficient. 

Solving  for  Z,  one  obtains  the  formula 

Z  »  LN(c,  ■  RJ/2kt  -  Xj/Zfcj, 

where  we  have  adopted  the  convention  that  X,  =  LN(L — 
L,„).  This  single-band  reflectance  model  assumes  that  the 
bottom  reflectance  is  constant  over  the  bottom  type,  that  the 
atmosphere  and  the  sea  state  are  uniform,  and  that  other 
background  optical  effects  are  either  uniform  or  constant 
throughout  the  image. 

To  reduce  errors  due  to  the  variation  of  the  bottom  reflec¬ 
tances,  a  two-band  ratio  method  (or  dual-band  method)  was 
devised.2’3  In  this  algorithm,  the  depth  is  given  by  the 
equation 

,  -M 

•  [LMc,  •  fi0l/c2  •  Ro2)  +  X,  -  X2], 


where  the  subscripts  1  and  2  indicate  different  bands.  In 
this  method,  the  assumption  is  that  changes  in  the  bottom 
reflectances  occur  so  that  the  ratio  Ci  •  /?„ i/c2  •  R„ 2  remains 
constant. 

Another  method,  which  we  call  the  linear  multiband  meth¬ 
od,  gives  the  depth  by 

where  the  sum  is  taken  over  several  bands  and  the  weights  ui, 
satisfy  the  constraint  that  2a>,  -  1.  Paredes  and  Spero,4 
generalizing  the  assumption  that  the  ratio  ci  •  R„t/c2  •  Ra2 
remains  constant,  assume  there  are  constants  ft  and  a,  inde¬ 
pendent  of  bottom  type  a,  so  that 

(ci  ’  Ra lH*l  ’  (c2  ‘  RaiR>2  '  (c3  *  ^03)1*3-  •  ■  =  a* 

With  this  assumption,  they  show  that  certain  weights  o>,  can 
be  found  so  that,  when  used  in  the  above  multiband  linear 
depth  equation,  produce  an  equation  for  Z  independent  of 
the  bottom  reflectance  and  depend  only  on  the  values  ft: 

Z  =  (V*  •  Zft  •  k,)  •  (1  -  f,  •  X,  -  fc  •  X2. . . ). 
Alternatively,  this  equation  may  be  written 

Z  =  A0+A1.X,  +  A2.X2+  ...  +  An  ■  Xn 

where  the  coefficients  A0,  At,  . . A„  are  constants  indepen¬ 
dent  of  the  bottom  type  at  which  the  depth  is  being  calculat¬ 
ed  (see  Ref.  4  for  more  details). 

It  has  been  suggested  that  when  the  assumption  underly¬ 
ing  the  two-band  ratio  method  is  not  satisfied,  the  linear 
multiband  algorithm  may  be  more  accurate  given  more 
bands  than  bottom  classes.  This  method  has  been  applied 
to  two  channels  of  multispectral  data  using  an  airborne  sen¬ 
sor.5  In  this  paper  we  use  the  linear  multiband  algorithm  on 
two  channels  of  a  Landsat  4  TM  scene  and  demonstrate  its 
improved  performance  over  the  two-band  ratio  method. 

The  Landsat  scene  (scene  5032614162  taken  21  Jan.  1985) 
contains  Isla  de  Vieques,  an  island  off  the  southeast  coast  of 
Puerto  Rico  in  the  Carribean  Sea.  The  water  clarity  in  this 
area  is  extremely  good,  and  bottom  reflectances  in  the  blue 
and  green  channels  (TM  bands  one  and  two)  are  detectable 
for  depths  up  to  25  m.  Bottom  reflectances  in  the  red 
channel  (TM  band  three)  are  detectable  for  depths  up  to  ~6 
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Fig.  1.  Depth  residuals  resulting  from  the  linear  multiband  model. 
The  residuals  were  calculated  using  the  second  set  of  calibration 
points  which  were  not  used  in  the  original  regression.  A  Gaussian 
distribution  fit  to  this  histogram  yielded  a  mean  of —0.13  m  and  a  a  ~ 
1.86  m. 


m;  sensor  noise  (striping)  over  the  water  renders  this  band 
useless  at  deeper  depths.  . 

The  general  procedure  for  this  comparison  was  to  first 
convert  a  NASA/GSFC  CCT  to  an  image  data  file  and  to 
georeference  this  data  file  against  a  Defense  Mapping  Agen¬ 
cy  1:25000  combat  chart.  This  combat  chart  contains  not 
only  bathymetric  soundings  but  also  land  feature  identifica¬ 
tions  that  aid  in  the  georeferencing  process.  A  set  of  ~600 
calibration  points  was  then  selected  by  recording  depths 
from  the  chart;  about  300  of  the  points  were  used  in  our 
regression  fit,  and  the  remaining  300  used  as  a  test  set  to 
check  the  calculated  depth  against  the  actual  depth. 

The  average  deep  water  signal  for  each  water  penetrating 
band  (in  this  case  channels  1  and  2,  roughly  corresponding  to 
the  blue-green  and  green  portions  of  the  visible  spectrum, 
respectively)  was  calculated  to  obtain  the  L,„  values.  These 
values  are  subtracted  from  the  corresponding  L,  value  to 
adjust  the  signal  values  for  atmospheric  scattering  etc. 

A  linear  regression  of  the  model  equation  against  the  (ap¬ 
propriate)  calibration  points  is  run,  which  produces  coeffi¬ 
cients  and  hence  equations  for  the  depth  Z  at  any  pixel.  The 
test  set  of  calibration  points  is  used  to  test  the  fit  of  this 
regression.  These  equations  are  then  used  to  produce  a 
bathymetric  image  that  can  be  processed  further  (smoothed, 
contoured,  pseudocolored,  etc.)  depending  on  the  desired 
application. 

The  linear  regression  was  run  on  the  equation 
Z  =  A0  +  AtXt  +  A^X^ 

corresponding  to  a  two-channel  multiband  linear  model. 
The  regression  yielded  the  following  values  for  the  constants: 
An  =  11.2;  A\  =  3.14;  A2  =  -6.76.  The  multiple  correlation 
coefficient  was  0.86  for  this  fit. 

Using  the  test  set  of  calibration  point  data,  the  model 
yielded  an  overall  residual  mean  of  <0.2  m  and  an  overall  rms 
of  <  1 .9  m.  These  residuals  are  shown  in  a  histogram  in  Fig. 

1.  In  the  depth  range  of  0-5  m,  the  rms  error  was  0.88  m;  in 


Fig.  2.  Depth  residuals  resulting  from  the  two-band  ratio  algo¬ 
rithm  calculated  using  the  second  set  of  calibration  points.  The 
residual  mean  is  0.093  m,  and  the  rms  is  2.66  m. 


the  5-10-m  range,  the  rms  error  was  1.25  m;  and  in  the  10-16- 
m  range,  the  rms  error  was  1.00  m. 

For  comparison  purposes,  the  dual-band  ratio  method  was 
applied  as  well.  To  minimize  error,  the  image  was  clustered 
to  locate  areas  of  similar  bottom  reflectance  using  a  super¬ 
vised  statistical  clustering  routine  and  the  maximum  likeli¬ 
hood  classifier.  The  algorithm  was  then  regressed  against 
calibration  points  in  each  cluster  separately.  An  overall 
residual  mean  of  0.093  m  and  an  rms  error  of  2.66  m  was 
obtained.  A  histogram  of  the  residuals  is  shown  in  Fig.  2.  In 
the  depth  range  of  0-5  m,  the  rms  error  was  0.95  m;  in  the  5- 
10-m  range,  the  rms  error  was  1.23  m;  and  in  the  10-16-m 
range,  the  rms  error  was  1.79  m. 

In  comparing  the  two  methods,  both  algorithms  underesti¬ 
mated  the  depth  in  shallow  water  and  overestimated  in  deep¬ 
er  water.  This  tendency  produces  a  larger  overall  error  than 
is  obtained  when  considering  the  indicated  depth  ranges.  It 
can  be  seen  that  the  linear  multiband  method  yields  some¬ 
what  improved  results,  even  when  only  two  bands  are  avail¬ 
able.  Moreover,  this  method  does  not  require  the  clustering 
and  classification  routines  to  discriminate  areas  of  similar  \ 
bottom  reflectance,  and  a  considerable  savings  in  CPU  pro¬ 
cessing  time  (one  to  several  hours  on  a  VAX  11/780  per  512- 
X  512-pixel  image)  is  realized.  * 

This  work  was  performed  at  the  Naval  Ocean  Research 
and  Development  Activity  and  sponsored  by  the  Office  of 
Naval  Technology,  Program  Element  62435N. 
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